Phase diagram of quantum capacitive junctions. 
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Mesoscopic capacitances, in the quantum regime, are described as a quantum rotor coupled to 
a dissipative bath. We analyze the phase diagram of this model, as function of the capacitance, 
the coupling to the bath (that is, the conductance of the device), and the decay of the correlation 
functions in the bath, i. e., the strength of the many body effects within the leads. By combining 
Monte-Carlo and analytical techniques, we present the full phase diagram. The similarities and 
differences with other models of quantum junctions are discussed. 

PACS numbers: 75.10.Jm, 75.10.Lp, 75.30. Ds. 
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Quantum dissipative systems Jjj3| are models for 
mesoscopic electronic devices. They consist of one, or 
more, tunnel junctions between normal or superconduct- 
ing electrodes. Quantum effects are associated to the dy- 
namics of a collective variable which describes the overall 
phase of the electronic wavefunction. For the case of su- 
perconducting systems, this variable is the phase of the 
order parameter |Q. A collective variable, conjugate to 
the total charge, can be defined for normal junctions as 
well [||J7| . The magnitude of the quantum fluctuations is 
governed by the capacitances in the device. 

Standard models describe these systems in terms of a 
few collective variables. In the simplest case, these vari- 
ables are the phase mentioned previously, and its conju- 
gate variable, the charge. The influence of the remaining, 
microscopic, degrees of freedom is taken into account by 
means of effective transport coefficients (viscosity, con- 
ductance) which couple the collective variable to a reser- 
voir, described in terms of harmonic oscillators. Two 
generic models have been analyzed within this frame- 
work: i) The Caldeira-Leggett model [Q which describes 
the dynamics of a junction where the quantization of 
charge is unimportant, and ii) The quantum rotor ||[|, 
in which only quantized charge transfer processes take 
place. The first model has been extensively studied. A 
generalization of it, the Schmid mod el p| , also describes 
junctions between Luttinger liquids [|10|. 

In the following, we will study in detail the second 
of these models. We consider a system made of gapless 
components, in which quantum (charging) effects are im- 
portant. We do not include, in the description of the 
environment, processes which destroy the charge quanti- 
zation in the device, such as an external electromagnetic 
field. We consider the effective action: 
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The model has four parameters: i) a short time cut 



off, t c , ii) the charging energy, Eq — (e 2 /C)~ 1 , where C 
is the capacitance of the device, iii) a dimensionless cou- 
pling, a ~ %/ (e 2 i?), where R is the resistance, and iv) the 
exponent g. The case of g = 1 corresponds to a quantum 
capacitance between two normal metals, j ^ 1 implies 
the existence of strong many body effects in the elec- 
trodes. They may arise because of shake-up and exciton 
effects induced by electrons as they tunnel, closely re- 
lated to the Mahan-Nozieres-de Dominicis (l^] processes 
in X-ray photoemission. An alternative way of inducing 
a value of g 7^ 1 is when the electrodes are one dimen- 
sional Luttinger liquids, a situation which can be realized 
experimentally |l3f| . 

The model given in ([!]), has been studied extensively 
by a variety of methods. A standard perturbative renor- 
malization group analysis can be performed when a>l 
and g ~ 1 jL4j. The scaling shows the existence of a 
phase transition along a critical line, a c cx (1 — g c ) ■ 

The case g — 1 has also been analyzed analytically 
jl5| [L7| ) and numerically |p^| , ^9| . The system shows a 
crossover from a high a regime, dominated by small am- 
plitude oscillations (spin waves), to a disordered phase 
for small values of a. The possibility of a sharp tran- 
sition, instead of a crossover, has also been postulated 
|]l6| , ^8 19|. Because of the periodicity in the variable 9, 
different sectors, characterized by the conjugate variable 
to 9, can be defined. The edge of the Brillouin Zone gen- 
erated in this way can be mapped onto a dissipative two 
level system [p|pc|], and solved. When Eq = 0, the rotor 
loses its dynamics. In terms of the environment, implicit 
in (|l|), we can define an effective Schmid model, with a 
non linear potential acting on the position closest to the 
rotor (see below). 

The model is also related to the "quantum rotor" prob- 
lem, introduced in the study of monopole induced proton 
decay in grand unified field theories J2IJ . This model can 
be solved exactly, in the limit Eqt c — » 0, (see |2l[] and 
also |2^]). As will be discussed later, this model is equiv- 
alent to eq. (|l|) for a — N and g — 1, where N is the 
number of fermion flavors coupled to the rotor. 
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We will study the action in eq. (|lj) by means of a nu- 
merical Monte-Carlo procedure, supplemented by ana- 
lytical methods, and by comparing to the extensive liter- 
ature mentioned above. The dynamics of the collective 
variable 8, arises from its interaction with the microscopic 
degrees of freedom of the environment. Instead of study- 
ing numerically a ID system with long range interactions, 
we analyze a 2D model with short range interactions. 
The environment is modeled by a 2D X-Y system, with 
couplings J x = J T = Jb, were the subscripts x and r 
stand for the space and the imaginary time dimensions. 
The planar spins are defined at the nodes of a square 
lattice, with periodic boundary conditions. We assign to 
each spin an angle, 4> X . T . In order for the environment 
to show algebraic correlations, we have Jb > J c ~ 1-12. 
For Jb 3> J c , a standard spin wave analysis gives: 

(cos(0 XjT - 4>x, T ')) oc 7c (2) 

(r - t') 4Jb 
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FIG. 1. Sketch of the model analyzed numerically in the 
text. 



in our discrete model constrains the effective interactions 
between the 9 T 's to be also periodic. The existence of al- 
gebraic decay in the environment (||), determines a sim- 
ilar decay for the effective retarded interactions of the 
rotor. An exact mapping of the parameters Jb,Jr,Jd 
into r Cl a,g can only be done numerically. A reasonable 
approximation can be obtained by integrating out the en- 
vironment using second order perturbation theory, and 
reexponcntiating the result. Then, as mentioned earlier, 
J 2 D . Finally, using the spin wave result (§), wc hnd 
2g = {n 2 J c )/{AJ B )- Note that, as J B > J c ~ 1.12, the 
only way of achieving values of 2g > 0.3 is to make n > 1 
(we will mostly take n = 4). 
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The rotor is defined as an additional ID X-Y model 
which is coupled to the environment in the way shown in 
fig. (0). It introduces two parameters, Jr and Ju- The 
first one determines the correlation length of the rotor 
when decoupled of the bath, and the second gives the 
magnitude of the interaction with the environment. In 
terms of the action (^), Jr oc -B^ 1 , and, approximately, 
Jjj) oc a. The rotor is described in terms of the angle 9. 
We couple it to the nearest column of environment spins 
by a term: 

Scoupling = Jd ^ C0s(6* t — n(f) XQ , T ) (3) 
r 

where n is an integer defined to facilitate the tuning of 
the parameter g in (^). 

By integrating out the environment degrees of freedom, 
we obtain an effective action, which, at long imaginary 
times, behaves like (|J). The periodicity of the variables 



FIG. 2. Phase diagram of the model. Circles are points 
analyzed by Monte-Carlo (Ec = 0.4, in umits of the cutoff). 
The full circles denote cases when an unambiguous character- 
ization of the correlation functions can be made. The lines 
denote the prediction of the critical line from the variational 
spin wave analysis discussed in the text. Full line, Ec = 0.4. 
Broken line, Ec = 0.2. Dashed line, Ec = 0.05. Arrows 
in the top part give the flow of the RG scaling equations for 
a S> 1 and g ~ 1. The two phases that we find are denoted 
C (capacitive) and D (dissipative) . 

We have studied the model by Monte-Carlo p3] . Typ- 
ical system sizes are l x = 20 — 40 and l T = 100 — 200. The 
number of configurations used is ~ 20000. The samples, 
produced by the Metropolis algorithm, were not corre- 
lated (autocorrelation < 0.1). The cases analyzed are 
shown in fig. (0). The two phases depicted in this figure 
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are characterized by the rotor correlation functions: 

G(t - t') = <cos(0 T - 6 r ,)) (4) 

This function shows two distinct regimes, as shown in 
fig. (H). For sufficiently large values of g, or low values of 
a, the correlation (|j) decays exponentially, while, if g is 
small, it shows an algebraic dependence on t - t'. 

The thick line shown in fig. (||) separates the two 
regimes. In terms of the original model, an algebraic de- 
cay in (|j) implies dissipative behavior, although it needs 
not be ohmic. When (|4|) decays exponentially, the sys- 
tem behaves as a capacitor, with an effective capacitance 
proportional to the inverse correlation length p3|. 



Semilogarithmic plot of the correlation functions 





g=l/4 - 




g=3/4 - 







distance n 

Logarithmic plot of the correlation functions 
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FIG. 3. Values of (cos(r — r')} as function of r — r'. 
The two cases shown correspond to g = 1/4 and a ~ 0.27 
(Jd/Jc = 0.5), full line, and g — 3/4 and a ~ 0.27, broken 
line. 



We now study the phase diagram of the model by vari- 
ational methods. The standard way is to employ the 
self consistent harmonic approximation (SCHA), which 
amounts to a spin wave approximation. The couplings of 
the model are replaced by quadratic interactions, whose 



strength is then optimized. This scheme has already been 
applied to the action (|l|) , for g — 1 jlj] , and at finite tem- 
peratures. For the related Schmid model, it reproduces 
correctly the phase diagram [ pif . 

For that purpose, we use the hamiltonian formalism: 



H = E c ^- + V [1 - cos (6- 0o)] + 



^ 2 

n=0 



(5) 



where [0,P] — [(f> n ,Pn] — i and a oc (V/uj c ) 2 . For 
a decoupled chain (V = 0), we have lim t ^ oo ([</> (t) — 
MO)?) = l/(47ro s )log(a; c *), so that 2g = 1/(4ttw s ). 
Note that, in this representation, the connexion of this 
model and the Schmid model in the limit Eq — > is 
transparent. 

The SCHA hamiltonian is: 



P 2 

T~Ischa — Ec — + E 



K ' 



-0o) 2 
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where TL c hain is the last part in eq. (||) and Ek 
is the variational parameter. We need to minimize 
(SCHA\H\SCHA), where \SCHA) is the ground state of 
TiscHA- hi the limit Ex, Eq *C uj c , we have: 



(SCHA\H\SCHA) = c^EcEk 



(7) 



where c is a numerical constant of order unity. 

Eq. (|7]) admits solutions for g < 1 and, roughly, 
a > c/(l — g)(Ec/i-iJcY 2 ~ 29 \ where c is a numerical con- 
stant of order unity. The factor (1 — g)^ 1 has been in- 
serted in order to interpolate correctly to the a — ► oo, 
g — * 1 limit. This line is plotted in fig.(||). 

The possibility of a description of the physics of the 
model in terms of only spin waves does not necessarily 
imply that other effects are irrelevant. Experience with 
related models, however, like the Schmid model ^Jj, or 
the ID quantum sine-Gordon gives some confidence on 
the reliability of the approach. 

We now turn to the relation of our work to the prob- 
lem of monopole induced baryon decay. The main issue 
addrssed in this context has been the dynamics of the 
baryons which build up the environment. It is, however, 
possible, to derive an effective action for the rotor coordi- 
nate. The leading terms should be those given in([l]). As 
the rotor is coupled to free fermions, g = 1. The analysis 
in the literature (2^j2^| , initially neglects the periodicity 
in the rotor coordinate. It implies, in expression (Q), the 
substitution of l-cos(0 T -0 T /) by (0 T - T ,) 2 /2. In Q, 
the coupling (that is, the fermion charge) is included in 
the definition of Ec, and the value of a is simply N, the 
number of fermion channels coupled to the rotor. This 
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model reduces to the Caldeira-Leggett model of a quan- 
tum particle in a dissipative environment Q), which can 
be solved exactly. The rotor coordinate experiences a 
frictional force, in agreement with the results in 

The solutions obtained within this scheme describe the 
model in terms of low energy electron hole pairs (spin 
waves, in our language). Such a description exhaust the 
physics of the problem in the right side of the phase dia- 
gram in fig. (H) ■ Thus, the validity of the scheme depends 
on the value of Ec/w c - It is assumed in [ plp^ ] that 
Ec /u> c — > 0, so that a spin wave description is marginally 
correct for g = 1. 

Note that a wavefunction which restores the symmetry 
of the model can be written by superposing \SCHA) and 
the solutions of H n — Hscha + It^uEk^ — 4>o)- Such 
state, e lkn \SCHA n ), has the same physical proper- 
ties as \SHA), because (SCHA n \SCHA n >) = 

It is interesting to speculate what would happen if 
E c /lo c is finite. Then, quantum fluctuations in the ro- 
tor coordinate would decouple it, at low energies, from 
the fermions. The rate of monopole induced baryon de- 
cay would not be of order unity, but would tend to zero 

Summarizing, we have studied the phase diagram of 
a model used to describe a variety of physical systems, 
such as Coulomb blockade in a normal metal junction, 
the single electron transistor or the quantum box. The 
generalization that we make in letting the parameter g in 
eq. ([j]) be different from 1 allows us to study other situ- 
ations, like charging effects in Luttinger liquid junctions, 
or shake-up and excitonic processes in normal junctions. 

We find a phase diagram with two regions: one in 
which charging effects dominate and the system behaves 
as a capacitor at low energies, and another in which the 
quantization of the charge is irrelevant, and the device 
shows a dissipative response. This phase diagram has 
been obtained by Monte-Carlo methods, and its consis- 
tency checked against alternative analytical approaches, 
such as the RG for a ^> 1 and g ~ 1 and variational 
ansatze for a -C 1. The phase diagram shows interest- 
ing deviations from the well understood phase diagram 
of the Schmid model: i) there is no duality, and the sepa- 
ratrix is not vertical, ii) there is a dynamically generated 
scale below which gaplike features appear (in the insu- 
lating phase) and iii) the charge quantization explicitly 
imposed in eq. (Q) reduces the the range of parameters 
for which the system shows dissipative behavior. 

We have studied the ground state only, and we have not 
considered the influence of finite voltages, which are im- 
portant in some situations, like the single electron tran- 
sistor. We can anticipate, however, that gate or bias 
voltages will only be relevant for the phase shown to the 
right in fig. (§). 

One of us (T. S.) thanks the Deutscher Akademischer 
Austauschdienst for financial support. This work has 



also been supported by the CICyT (Spain) through grant 
MAT94-0982. 



[1] 
[2] 

[3] 
[4] 
[5] 
[6] 
[7] 



[9] 
[10] 

[11] 
[12] 



[13] 

[14] 
[15] 
[16] 

[17] 



[18] 
[19] 

[20] 

[21] 
[22] 



A. O. Caldeira and A. J. Leggett, Ann. of Phys. (NY) 
149, 374 (1983). 

D. V. Averin and K. K. Likharev in Mesoscopic Phe- 
nomena in Solids, edited by B. Altschuler, P. A. Lee 
and R. A. Webb (Elsevier, Amsterdam, 1991). 

U. Weiss, Quantum Dissipative Systems, World Sci- 
entific, Singapore (1993). 

V. Ambegaokar, U. Eckern and G. Schon, Phys. Rev. 
Lett. 48, 1745 (1982). 

E. Ben- Jacob, E. Mottola and G. Schon, Phys. Rev. Lett. 
51, 2064 (1983). 

F. Guinea and G. Schon, Europhys. Lett. 1, 585 (1986); 
J. Low Temp. Phys. 69, 219 (1987). 

G. Schon and A. D. Zaikin, Phys. Rep. 198, 238 (1990). 
G. Schon, in Quantum Processes and Dissipation, 
Springer (Berlin), in press. 

A. Schmid, Phys. Rev. Lett. 51, 1506 (1983). 

C. L. Kane and M. P. A. Fisher, Phys. Rev. Lett. 68 1220 

(1992). Phys. Rev. B 46, 15233 (1992). 

M. Ueda and F. Guinea, Z. Phys. B 85, 413 (1991). 

G. D. Mahan, Many-Particle Physics, Plenum (New 

York) 1991. P. Nozieres and C. T. de Dominicis, Phys. 

Rev. 178, 1097 (1969). 

F. P. Mulliken, C. P. Umbach and R. A. Webb, as re- 
ported in Physics Today 47, 19 (1994). 

J. M. Kosterlitz, Phys. Rev. Lett. 37, 1577 (1977). 
R. Brown and E. Simanek, Phys. Rev. B 34, 2957 (1986). 
W. Zwerger and M. Scharpf, Zeits. fur Phys. 85, 421 
(1991). 

S. V. Panyukov and A. D. Zaikin, Phys. Rev. Lett. 67, 
3168 (1991). D. S. Gobulev and A. D. Zaikin, Phys. Rev. 
B 50, 8736 (1994). 

R. Brown and E. Simanek, Phys. Rev. B 38, 9264, (1988). 
V. Scalia, G. Falci, R. Fazio and G. Giaquinta, Phys. 
Rev. Lett. 67, 2203 (1991). 

G. Falci, G. Schon and G. T. Zimanyi, Phys. Rev. Lett. 
74, 3257 (1995). 

J. Polchinski, Nuclear Physics B242, 345 (1984). 

I. Affleck and J. Sagi, Nucl. Phys. B417, 374 (19 94). J 

M. Maldacen a and A. W. W. Ludwig, preprint (cond 

mat/9502109|). 



[23] The spectrum of the system, however, is gapless, as there 
is no gap in the environment. The low energy modes of 
reservoir always couple to the rotor, influencing its long 
time dynamics. See F. Guinea and M. Ueda, Phys. Rev. 
B 49, 5722 (1994). 

[24] M. P. A. Fisher and W. Zwerger, Phys. Rev. B 32, 6190 
(1985). 

[25] K. Binder, D. W. Heermann, The Monte-Carlo 
Method in Statistical Physics, Springer Series in 
Solid State Sciences, 80, Springer (Berlin, Heidelberg), 
1988. 



4 



